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AN ANALYTICAL EVALUATION OF THE DENTING 
OF AIRPLANE SURFACES BY HAIL 

By Robert G. Thomson and Robert J. Hayduk 
Langley Research Center 

SUMMARY 

Hailstone damage to airplanes is exhibited as denting or complete penetration. An 
analytical evaluation of such denting characteristics of airplane surfaces is idealized by 
having a crushable spherical projectile impact normal to a flat plate. Separate elastic 
and rigid-plastic analyses have been developed previously by using the classical thin- 
plate bending approach. As applied in this paper, the elastic analysis determines the 
area of the sheet around the center of impact, which has become plastic. The rigid- 
plastic analysis applied to this central plastic region determines the resulting permanent 
deformation (dent). Numerical results are presented in detail for the permanent defor- 
mations of two specific sheet thicknesses. 

INTRODUCTION 

Because of the high velocity of modern airplanes, the probability of surface 
damage caused by impact with hail, and possibly with other meteorological precipitants, 
has become a design consideration. Meteorological damage can take the form of either 
surface -material erosion, permanent deformation, or, if the impact is severe enough, 
complete penetration of the surface. A summary of the effect of hail on airplanes in 
flight with photographic documentation of damage and an extensive bibliography is given 
in reference 1. 

No design criterion for hail encounter has been available to the aircraft designers 
until recently. In reference 2 a testing procedure is proposed to aid the designer in 
minimizing hail damage. The suggested approach is to test the affected structure by 
impacting it with representative samples of hail at various impact velocities to assure 
that no appreciable damage occurs. An analytical approach does not appear to be avail- 
able, and the purpose of the present paper is to present an approximate analytical method 
to complement established laboratory testing procedures. 

When a hailstone impacts the thin skin of an airplane, the resulting skin deforma- 
tions usually are plastic near the point of impact and elastic in the surrounding area. If 



the impact is severe enough, both bending and membrane action are present in the skin. 

A rigorous, but difficult, approach to the problem would require a sophisticated, dynamic, 
elastic -plastic, large-deformation analysis. The analytical development in this paper, 
however, attempts to approximate the analysis of the problem to circumvent the difficul- 
ties of the more rigorous approach. 

In reference 3 an elastic analysis was developed by using the classical thin-plate 
bending approach for the impact of a spherical projectile on a thin infinite sheet. The 
results of that analysis determine the time -dependent expressions for the deformation, 
velocity, and moment distributions of the elastically deforming sheet. The elasticity 
assumptions, however, preclude any permanent -type deformations. 

The rigid-plastic bending response of a finite circular sheet subjected to a trans- 
verse impulse loading has also been investigated. (See ref. 4.) The rigid-plastic analy- 
sis considers the entire finite sheet as plastic, and when deformations occur, the whole 
sheet experiences a permanent type of deformation. Therefore, in regard to hailstone 
impact, if the extent of the plastic area around the center of impact is known, a rigid- 
plastic analysis can be used to predict the permanent-deformation characteristics 
(denting) of this region. 

In the present paper the elastic analysis of reference 3 is used to determine the 
radial and circumferential bending-moment distributions in an infinite elastic sheet 
impacted by a crushable spherical projectile. A rigid-plastic analysis is also developed 
herein for a simply supported finite circular sheet impacted by a crushable spherical pro- 
jectile, and the resulting permanent deformation is found. The two analyses have identical 
initial conditions on loading: a crushable spherical projectile in which the projectile is 
assumed to impose an initial, instantaneous, finite velocity distribution over the region of 
impact, the rest of the plate being at rest. In both analyses the circumferential and radial 
moments are at a maximum at the center of the impact and decrease monotonically with 
increasing radius. For the infinite elastic sheet the bending moments around the center 
of impact are shown to exceed the dynamic elastic yield-moment resultant M 0 . 

In the present approach the simply supported, rigid-plastic boundary conditions on 
the radial and circumferential bending moments of the plastic analysis are matched with 
those of the elastic analysis at some radius near the center of impact called the plastic 
radius. The rigid-plastic analysis, using this plastic radius as the radius of the finite 
circular sheet, is then to be used to determine the permanent deformations. The principal 
contention of the present paper is that sufficiently accurate results may be obtained when 
only the boundary conditions on the bending moments are satisfied between the elastic and 
plastic analyses in determining the plastic radius.' 
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The permanent deformations on the rigid-plastic sheet increase with increasing 
hailstone size and velocity. However, it is recognized that a failure by perforation occurs 
at a certain critical impact velocity and that this failure results in a practical limit to the 
deformations. This practical limit can be established once experimental data become 
available on the dynamic material properties associated with hailstone impact. A pre- 
liminary estimate of the location of this limit is presented in this paper. 

SYMBOLS 

The units used for the physical quantities defined in this paper are given both in U.S. 
Customary Units and in the International System of Units (SI). (See ref. 5.) Appendix A 
presents factors relating these two systems of units. 

a radius of projectile 

b radius of central plastic region 

c = j/E/p^ 

D flexural rigidity of sheet, Eh^/l2 (l - vty 

E elastic modulus of sheet material in tension and compression 

g Q projectile velocity 

h target thickness 

J n (z) Bessel function of the first kind 

k r ,kg radial and circumferential time rates of curvature 

M bending-moment resultant 

M 0 yield-moment resultant, cr 0 ( h/2)^; maximum bending moment that sheet 

can sustain 

m mass 
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Hankel transform parameter 


transverse shear -stress resultant 


impulse loading 


radial, circumferential, and transverse coordinates, respectively (see fig. 1) 


time of cessation of all sheet motion after impact 


time at which hinge -circle radius has decreased to zero 


2a 

h P t 


sheet deflection in transverse direction (z-direction) 


permanent center deflection, w(0,t*) 


radius of hinge circle separating regimes A and AB 


mass per unit area of sheet, pj.h 


Poisson's ratio 


mass density 

yield stress in simple tension or compression (assumed to be identical in 
magnitude) 


V y ,Oq 


principal radial and circumferential stresses, respectively 
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r = t/k 


i// momentum constant (see eq. (7)) 

Subscripts: 
cr critical 

n indicates order of Bessel function 

p projectile 

r radial 

t target 

9 circumferential 

A dot over a symbol indicates partial differentiation with respect to time t. An 
asterisk denotes time at which the entire sheet comes to rest. 

MATHEMATICAL ANALYSES 

Determination of Initial Velocity Distribution 

Linear, small-deflection plate theory was used in reference 3 to study the response 
of an infinite elastic sheet subjected to an axisymmetric impulse. The analysis assumed 
that the momentum exchange was the primary loading mechanism, that the exchange was 
instantaneous, and that the mass of the projectile was negligible after impact compared 
with the mass of the sheet. The three specific projectile shapes studied in reference 3 
were cone, sphere, and cylinder. Reference will be made in this paper to the solution 
for the spherical projectile only, and the momentum distribution shown in figure 1 is 
derived from this shape. 

The assumed instantaneous momentum exchange between projectile and sheet yields 
an initial sheet velocity over the centrally impacted area, but no initial displacement. In 
the present analysis the mass of the projectile after impact is not neglected as in refer- 
ence 3. The momentum balance for a completely inelastic impact is 


( 1 ) 


(m p + m t )w( 7 j, 0 ) = m p g Q 
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or, from figure 1 


27ra 2 r?(2a^l - r? 2 )p p dr? + 27ra 2 r?hp t d?7 w(r?,0) = 2Tra. 2 ri(pa.)Ji - r/ 2 )p p dr? 


So 


( 2 ) 


Solution of equation (2) for the initial velocity distribution w(r?,0) gives 


w(?7,0) 


2a 

h 


1 + 



w(?7,0) = 0 


(0 ^ 77 £ 1) 

> (3) 

(V = 1 ) 


An approximation to this initial velocity distribution is made to utilize the results 
of reference 3 in which the initial velocity is assumed to be of the form: 



w(?7,0) = 0 


(0 ^ 77 = 1) 

(r? £ 1) 


(4) 


In this approximation the character of equation (4) is retained but the mass of the pro- 
jectile after impact is included. The approximation to equation (3) is 


w(r?,0) = V 0 {l - r? 2 
w(?7,0) = 0 


(0 s I? < 1) I 

(V = 1 


(5) 


where 


V o = 


Pt 


bu 




h p t 


The momentum constant \p is determined by equating the total momentum of the sheet 
for two sheet initial velocity distributions. In one case, the exact velocity distribution 
(eq. (3)) is employed in the momentum equation. In the other case, the approximate 
velocity distribution (eq. (5)) is used; that is, 
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( 6 ) 



/l - 7)2 

1 + IT 
h Pt 


rj d?7 


and the momentum constant i// becomes 


3. 3_ 3 ln 4 , 2aPp\ 2a Pp 

2 2a Pg /2a Pp\ 2 \ h p tJ h p t 

h p t \ h P t J 

The values of ^ used in the specific calculations in this paper are presented in table I 
along with the target (2024 -T4 aluminum) and projectile material (ice) properties. A 
comparison of the exact initial velocity distribution (eq. (3)) with the approximate initial 
velocity distribution (eq. (5)) is shown in figure 2 for a 0.25-inch-radius (0.64-cm) hail- 
stone impacting a 0.063 -inch-thick (0.1 60 -cm) plate (\p = 0.71). 


Determination of Moments in the Elastic Region 

The boundary conditions of zero deflection, slope, and shear per unit length at 
infinity and of zero slope and finite shear at the origin were assumed in reference 3 in 
order to obtain a solution in single integral form for the deflection of the dynamically 
deforming plate. The deflection (eq. (31) of ref. 3) can be written as 

W = k ir(?f g o)/l j* 0 P’ 5/2j 3/2(P) J 0(PP) sin(p2r)dp (8) 


2 2L 

where in the present paper ^“p^go will be replaced by V Q . The expressions for the 
radial and circumferential moments from elasticity theory are, respectively, 


M r (v,t) = ~ 


D 


9 2 w v 9w 

9p2 Tj dr] 


(9) 


M 0(77,7) = - 


D (l 9w 
a 2 V? dr l 



( 10 ) 


The radial and circumferential moments for v = 1/3 (eqs. (8), (9), and (10) combined) 
can be written in the form (see appendix B): 
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DkV rtl g 


M r (r 7 ,T) = a 2 2 J 3/2 (p) sin(p 2 r) p -1 / 2 |3 J 2 (pv) - Jq(p^) 


- |p _ 1 Jl(P»7)> dp (11) 


DkVj? 


M 0 (?7,t) 6a^ i() J3 / 2(P) Sin ( p2T ) J l(P^ + J 2(P?l) " JoCP 7 ?) 

and, for 77 = 0, they reduce to 

2DkV 0 i/| r oo 1/0 

M(0,r) = — — J 0 J 3/2(P) sin(p 2 T) p _1/2 dp 

M r (?7,T) 


dp (12) 


(13) 


The nondimensional radial and circumferential moment distributions 


M 0 (?7,t) 


M(0,t) 


and 


, respectively, in the central portion of the infinite elastic sheet (center of impact) 


M(0,t) 
are plotted in figure 3. 


Rigid-Plastic Analysis 

In reference 4 the dynamic rigid-plastic response of a thin, circular plate to an 
axisymmetric impulse loading with radial variation was considered. Specifically, the 
radial impulse loading had a Gaussian distribution. The circular sheet was assumed to 
be finite in extent, simply supported along its outer edge, and uniform in thickness. The 
sheet was also assumed to be composed of a material that is incompressible, ideally 
plastic (no work-hardening), and isotropic in yielding. Only the bending action of the 
plate was considered; the elastic deformations as well as the membrane forces were 
neglected. The rigid, perfectly plastic plate material was assumed to deform according 
to the Tresca yield condition and the associated flow rule. (See ref. 4.) 

The present analysis is similar to the analysis of reference 4. However, the axi- 
symmetric loading function is no longer Gaussian but is replaced by a projected spherical 
impulse distribution, and the effect of the mass of the projectile is not neglected after 
impact but is accounted for approximately as in the previously discussed elastic analysis. 
Hence, the rigid-plastic plate, which is initially at rest, immediately upon impact assumes 
a velocity distribution as given by equation (5); and both the elastic and rigid-plastic plate 
analyses have identical initial conditions: 

w(? 7 , 0 ) = 0 

w(i7,0) = V 0 /l - r? 2 (0? n SI) 

w(? 7 , 0 ) = 0 (77 = 1) 
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The details of the rigid-plastic solution are given in appendix C. The permanent defor- 
mation of the rigid-plastic plate can be written as 


w*(4,t*) 


pV 0 2 b 2 

24M 0 




24M q 

MV 0 b 2 


t(S) + (1 - & 


1 + 2£ + 3£ 2 




w*($,t*) 


M V 0 2 b 2 

24M 0 


(1 - I) 





(15) 


where the permanent center deformation 6 is expressed as (see fig. 1) 


5 = w*(0,t*) 


MVp 2 b 2 / 24 M q \ 

24M o ^ + M v o b 2 J 


(16) 


and £ = r/b, b is the radius of the plastic plate, and a is the radius of the projectile. 

24m 

The nondimensional time parameter — ■ -° t^ (see eq. (C9)) for 0 = a/b = 1 is 


pv 0 b 2 


24M n 

pv 0 b 2tl 


= -2 + 






(17) 


and for a/b ^ 1 (that is, projectile larger than plastic dent) is (see eq. (CIO)) 


24M 0 

MV 0 b 2 


tl 




(18) 


Also, the nondimensional time parameter 


24M 

%t(£) is expressed as (see eq. (Cll)): 

pV 0 b 2 


24M 0 

MV 0 b 2 


t(*) = 




- 2 £ + 2| 2 


/a\ 2 1 

- 3 (s) 




+ 


24M 0 


(19) 
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RESULTS AND DISCUSSION 


Determination of Radius of Plastic Center 

The simple -support boundary conditions of the rigid-plastic sheet, according to 
plasticity theory (see appendix C), are that the radial moment M r equal zero and the 
circumferential moment Mg equal M 0 . Hence, the particular radial location b/a of 
the central plastic portion of the infinite elastic sheet was taken to be at the radial loca- 
tion where M r = 0; this radial location changes with time. For different specific times 
the value of the circumferential moment Mg at the radial location where M r = 0 will 
possess values that range from above to below the yield moment M 0 (assuming that the 
sheet and hailstone properties remain constant). When Mg = M 0 (with M r = 0), and 
only then, are the boundary conditions met for the simply supported plastic sheet; and it 
follows that this radius is the only radius of plastic yielding on the elastic sheet at which 
the boundary conditions are satisfied. 

Typical (r = 1.5) radial (eq. (11)) and circumferential (eq. (12)) moment distribu- 
tions in the central portion of the infinite elastic sheet are plotted in figure 3. The maxi- 
mum (and only) radius of the rigid-plastic region for a particular hailstone-sheet combi- 
nation is determined inversely by matching the rigid-plastic boundary conditions at a 
particular radius of the infinite elastic sheet. For example, in figure 3 b/a is chosen 
to be at the radial location r) = 2.1 A where M r = 0; the value of Mg at rj = 2.74 
depends on the sheet and hailstone properties. If the sheet and hailstone properties are 
then selected so that Mg equals M 0 , the yield moment, then the radius of the rigid- 
plastic sheet b/a equals 2.74 for this particular hailstone -sheet combination. Hence, 
any elastic moment distribution specifies the radius of a rigid-plastic sheet for a particu- 
lar hailstone diameter, plate thickness, and material when the hailstone velocity is 
selected to make the elastic circumferential moment Mg (eq. (12)) equal to the yield 
moment M 0 . 

By calculating various moment distributions at a number of specific times, such as 
the one shown in figure 3, and applying the criterion for determining the maximum radius 
of plastic yielding (shown as the dashed line in fig. 3), a relationship between b/a and 
the initial projectile velocity can be established. Such a relationship is shown plotted in 
figure 4 for three different hailstones of 0.50-, 0.25-, and 0.125-inch radius (1.27-, 0.64-, 
and 0.32 -cm). In figure 4 two typical aircraft skin thicknesses were chosen as target- 
plate thicknesses: 0.04 inch (0.10 cm), shown as the solid line, and 0.063 inch (0.160 cm), 
shown as the dashed line, respectively. The nondimensional plastic radius is shown in 
figure 4 to increase, for a particular velocity, as either the target thickness decreases 
and/or the projectile radius increases. Also, the plastic radius increases as the projec- 
tile velocity increases, as is expected. 
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At extremely small values of t (that is, r 2= 0.2) corresponding to values of 
b/a > 1, it should be noted, however, that an anomaly occurs in the moment distributions. 
Since it is known that classical plate theory predicts a spurious response at all points of 
a plate immediately after the application of a sharp impulse (ref. 6), the anomaly in the 
moment distributions is felt to be mathematical rather than physical. The anomaly tends 
to distort the smooth monotonically decreasing form of the moment distributions as pre- 
sented in figure 3 and also the computed values of the plastic radius as shown in fig- 
ure 4. The computed values of the plastic radius begin to oscillate as b/a decreases 
below 0.5, and the curves presented in figure 4 have been faired between b/a = 0.5 and 
b/a = 0. The value of the projectile velocity at b/a = 0 can be easily (and accurately) 
determined by setting M(0,t) in equation (13) equal to M 0 . Fortunately, this oscilla- 
tion in the computed values of b/a occurs in the region where the plastic dent has a 
smaller radius than the impacting hailstone and generally is of less interest than in the 
region where b/a > 1. 

Determination of Permanent Center Deformation 

The radius b, which is the maximum radius of plastic yielding from the elastic 
analysis, is now to be used as the finite radius of the rigid-plastic plate. By use of equa- 
tion (16) and the appropriate values of a, b, h, M 0 , /i, and others (see table I) for 

target and projectile material properties, the permanent center deflection 6 of the 
plastically deformed plate is determined. In figure 5, a logarithmic plot of the permanent 
center deformation as a function of projectile velocity is shown for three different hail- 
stones of 0.50-, 0.25-, and 0.125-inch radius (1.27-, 0.64-, and 0.32-cm) and with a 
target-sheet thickness of 0.063 inch (0.160 cm). The permanent center deformation is 
seen to approach a linear asymptote (shown as a dashed line) as the projectile velocity 
increases for the three sizes of hailstones. This linear asymptote has a constant slope 
of 2.85. The permanent center deformation is also shown to be very sensitive to projec- 
tile size, probably because of the thinness of the target sheet relative to the mass of the 
hailstones. 

In figure 6 the shape of the permanent deformation is plotted (eqs. (14) and (15)) for 
three specific projectile impact velocities of the 0.50 -inch-radius (1.27-cm) hailstone. 
The projectile velocities are 330, 500, and 1100 ft/sec (0.10, 0.15, and 0.34 km/s) and 
show a progression from a shallow, circular, dish -shaped dent to a deep conical cavity. 
Thus, the shape of the dent, as well as the center deformation, changes as the projectile 
velocity is varied. The limited damage information available (ref. 1) indicates that the 
indentation predictions of the analysis are reasonable; however, adequate data for com- 
parison are not available. 
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In figure 7, a logarithmic plot of the permanent center deformation as a function of 
projectile velocity is shown for hailstones of 0.50-, 0.25-, and 0.125-inch radius (1.27-, 
0.64-, and 0.32-cm) but for a target-sheet thickness of 0.04 inch (0.10 cm). The perma- 
nent center deformation is again seen to approach a linear asymptote as the projectile 
velocity increases, in fact, the same asymptote as in figure 5. From the results given in 
figures 5 and 7 it is apparent that the permanent center deformation, in the high velocity 
range, is a function of the projectile velocity g Q raised to the 2.85 power. In an attempt 
to show the six curves given in figures 5 and 7 on a single curve, the permanent center 
deformation is plotted as a function of g 0 a/h2/3 in figure 8 to a logarithmic scale. The 
curves in figure 8 show a general tendency toward a common linear asymptote in the high 
velocity range, but some spread does exist among the curves. This spread can be attrib- 
uted, in part, to the choice of the abscissa parameter g 0 a/h^/^. 


Critical Impact Velocity 

The critical impact velocity is defined as the minimum impact velocity for a given 
sheet and projectile combination that barely allows the projectile to perforate the sheet. 

An analytical expression for the critical impact velocity of a thin sheet impacted by a 
spherical projectile is presented in reference 7. It was derived from the elastic response 
of the sheet by assuming a failure criterion based on a maximum fracture stress. The 
critical impact velocity of a thin visco-plastic sheet impacted by a cylindrical projectile 
was studied in reference 8. By comparing the results presented in references 7 to 10, it 
can be shown that the relation between the critical sheet thickness and the critical impact 
velocity is linear, as is reported in both the elastic and visco-plastic analyses. The fol- 
lowing analytical relationship (ref. 7) is felt to be a reasonable statement of proportion- 
ality to determine the critical sheet thickness h cr for thin sheets impacted by spherical 
projectiles: 


^cr 



( 20 ) 


A proportionality constant, which would be a function of the sheet material properties, can 
be determined once meaningful experimental data on the dynamic material properties of 
thin sheets impacted with hailstones become available. Limiting ballistic envelopes are 
plotted in figures 5 and 7 as the long-short, dashed curves. The shape of this limiting 
envelope is presumed to be known at present; its location has been estimated and is sub- 
ject to change. 
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CONCLUDING REMARKS 


The present analysis, although not an exact elastic -plastic treatment, has combined 
separate elastic and rigid-plastic analyses to study hail damage to airplane surfaces. The 
results of this analysis indicate that thin structural members, such as the large skin areas 
of airplanes are quite susceptible to denting by hailstones. The larger the hailstone the 
greater the damage, of course, but the analysis also indicates great sensitivity to impact 
velocity as the hailstone size increases. In fact, the permanent center of deformation of 
an impacted skin has been shown to approach a function of the impact velocity to the 
2.85 power for center deformations of 0.10 inch (0.25 cm) or over. 

The analytical results presented herein must await experimental verification. 
Although the magnitudes of the permanent deformations may be found to be in need of 
adjustment because of the neglect of membrane action and the numerous simplifying 
assumptions of the present analysis, it is hoped that the parametric study made herein 
and especially the impact-velocity dependency uncovered may prove to be useful. It is 
realized also that the results presented are for normal impact onto flat sheets, and that 
some of the most severe hail damage to aircraft occurs at the leading edge of wing and 
tail assemblies, which are highly curved. The strengthening and geometric factors con- 
nected with curvature of the sheet need to be evaluated. It is suggested that for oblique 
impact the normal component of the impact be considered as a first approximation. 

Perforation of the airplane skin occurs at some critical impact velocity, but the pre- 
diction by analysis of the magnitude of this impact velocity awaits experimental determina- 
tion of some material properties at impact conditions. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., May 15, 1969, 

126-61-01-07-23. 
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APPENDIX A 


Mill Hill 


CONVERSION OF U.S. CUSTOMARY UNITS TO SI UNITS 


The International System of Units (SI) was adopted by the Eleventh General Confer- 
ence on Weights and Measures, Paris, October 1960, in Resolution No. 12 (ref. 5). Con- 
version factors for the units used herein are given in the following table: 


Physical quantity 

U.S. customary 
unit 

Conversion 

factor 

(a) 

Length 

in. 

0.0254 

Density 

lbm/ft 3 

16.02 

Stress 

psi (lbf/in2) 

6.895 x 10 3 

Velocity 

ft/sec 

0.3048 


SI unit 


meters (m) 

kilograms per cubic meter (kg/m 3) 
newtons per square meter (N/m 3 ) 
meters per second (m/s) 


^Multiply value given in U.S. customary unit by conversion factor to obtain equiva- 


lent value in SI unit. 


Prefixes to indicate multiple of units are as follows: 


Prefix 

Multiple 

centi (c) 

10-2 

kilo (k) 

10 3 

giga (G) 

109 
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APPENDIX B 


MOMENT EXPRESSIONS FOR CRUSHABLE SPHERICAL PROJECTILE 


The expressions for the radial and circumferential moments of an elastic plate from 
elasticity theory are, respectively (see eqs. ( 9 ) and ( 10 )), 


,, In r\ D f a 2 W , V 3 w\ 




(Bl) 

(B 2 ) 


The elastic deflection in single -integral form is (eq. (8)) 


w = kv o/§ J 0 P _5//2 J 3/2(P) J 0 (P 7 ?) sin(p 2 r) dp 


(B 3 ) 


Hence, the radial and circumferential moments (by substitution of eq. (B 3 ) into eqs. (Bl) 
and (B 2 )) become 

M r (v,r) -gS. J 3 / 2 (P) sin(p 2 r)||p' 1 / 2 [j 2 (pr?) - Jo(P^)J- ^P _3,/2 Jl(P dp (B4) 

M e (v,r) = J3/2(P) sin(p 2 r)^2_Z_ j 1 ( p7? ) + |p-l/ 2 jj 2 ( p?7 ) - J 0 (prj)]J> dp (B 5 ) 

The evaluation of M r and at 77 = 0 is determined by noting that 


lim 
77 0 


lim 

77—O 


dJl(Pf?) 

drj 


Jq(P 7 ?) " J 2 (P r ?) 


'77 — O 


= P 
2 


Jl(P^) _ p2(P T ?) +Jq(PV) 


V 


1 77 — 0 


_ p 
2 


(B6) 


(B 7 ) 
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APPENDIX B 


If v is taken to be equal to 1/3, equations (B4) and (B5) reduce to 


M r (0,r) = M 0 (O,t) = M(0,t) = § J 3 /2(p) sin(p 2 r) p^/2 dp ( B8 ) 

^ 0 

The nondimensional radial and circumferential moments, respectively, can then be 
written as 


M r (77 


T ) So J 3/2^ sin (p 2r ) P" 1 / 2 j^Cpr?) ~ J 0 (Pn)] -§P _1 J l(P^)J d P 


M(0,t) 


J 3/2 (p) sin (p 2r ) P 1//2 d P 


(B9) 


M(0,t) 


"4 S 0 J 3/2(P) sin(p 2 r) p _ 1 / 2 ^-|p _1 J 1 (p??) + \j2(PV) ~ Jo(PT?)]} d P 

r OO 

J J3/2(p) sin (p 2r ) P " 1 / 2 dp 


(BIO) 
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APPENDIX C 


RIGID-PLASTIC SOLUTION FOR SPHERICAL PROJECTILE IMPACT 

The analytical derivation of the deflection of the plastically deforming, simply sup- 
ported circular plate presented in this appendix is similar in approach to existing dynamic 
rigid-plastic solutions. The development of the governing equations closely parallels that 
of reference 4, except for a change in the functional form of the initial velocity distribu- 
tion across the plate. (For a more complete discussion of rigid-plastic theory of plas- 
ticity, utilizing the Tresca yield condition and the associated flow rule, the reader is 
referred to reference 4.) The development of the deflection expressions for the plasti- 
cally deforming plate is divided in this appendix into two sections. The first section pre- 
sents the solution for the motion of the plate when the plate is divided into two distinct 
plasticity regimes, separated by a hinge circle, up to the time of the disappearance of 
this hinge circle. The second section contains the solution for the motion of the plate as 
it continues to deform, up to the time of cessation of all motion. Both sections will be as 
brief and concise as possible because the analytical development is similar to that pre- 
sented in reference 4, and a more complete treatment is not deemed to be necessary. 

Solution With Hinge Circle of Finite Radius 

The rigid-plastic plate, which is initially at rest, immediately upon impact assumes 
an axisymmetric velocity distribution in the central portion of the plastic plate (designated 
regime A of the Tresca yield hexagon in fig. 9) of the form 


w(77,0) 




(Cl) 


where the magnitude of the velocity V Q is a function of the projectile and plate param- 
eters (see eq. (5)), 


v 0 = 


2a P p 
h p t g o 

hPt 


and £(t) is the radius of the hinge circle separating regimes A and AB, and £ = r/b. 

In the surrounding annulus of the remaining portion of the plastic plate (designated 
regime AB of the Tresca yield hexagon in fig. 9), the proper velocity distribution 
according to the dictates of rigid-plastic theory (see ref. 4) is linear in £, continuous 
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across the plate, and zero at the simply supported edge. The satisfaction of these condi- 
tions leads to the following expression for the velocity distribution in regime AB: 


*<^.°> = v c f - (|) 2 (£) 2 (5 $ « = ») < C2 > 

1 “ b 

From the symmetry of the impulse loading it follows that the only nonvanishing stresses 
(on a circular plate element) are a r and oq and the transverse shear -stress result- 
ant Q r . Summation of the vertical forces and summation of the moments acting on the 
plate result in the following equation of motion: 


d_ 


(£M r ) = M 0 + b 2 y Oiw-q)£d£ 


where the impulse loading q is prescribed by 


(C3) 


q = M V 0 ^/l - 6(t) 

Equation (C3) is automatically satisfied by equation (Cl) and reduces to /l tw = q 
since M r = M q = M 0 in regime A of the Tresca yield hexagon (see ref. 4). In 
regime AB the equation of motion (eq. (C3)) is not uniquely satisfied until the radius of 
the hinge circle £(t) is determined. In order to solve for the proper time variation of 
£(t), the equation of motion is utilized with the velocity expression (eq. (C2)) and the asso- 
ciated boundary conditions for the surrounding annular region to obtain the following 
first-order total differential equation in ?(t)/b: 


*2M 0 

M V 0 b 2 


dt = 




1 - I bf li 2 

cl 


172 




4 


(C4) 


where the associated boundary conditions of regime AB are 


M r (|»t) = M 0 


M r (l,t) = 0 


(C5) 
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Integration of equation (C4) yields 


12M 0 

MV G b2 


t + Cj 



3/a 
2 lb 


£ 

b 



sin 


-1 £/b 
a/b 


(C6) 


where Cj is the constant of integration. 

The initial location of the hinge-circle radius £(0)/b is determined from two 
requirements: (1) that the time rates of curvature k r and kg be nonnegative in 

regime A, and (2) that d(£/b)/dt (eq. (C4)) be negative. The first condition of nonnega- 
tive time rates of curvature is satisfied by the initial velocity distribution in regime A for 
0 S | = a/b. (See eq. (Cl).) However, if the initial hinge circle radius is taken equal to 
a/b, then d(£/b)/dt would be zero, would indicate no movement of the radial position of 
?(t)/b with time, and no decrease in the total kinetic energy of the sheet and hence, no 
plastic deformation. If the position of ?(0)/b is taken to lie between (a/b)2 and a/b 
the rate of hinge-circle movement, d(£/b)/clt, is positive and the radial position of £(t)/b 
increases with time until it reaches the position a/b where the hinge-circle movement 
ceases once again. If, however, the maximum value of £/0)/b which yields a negative 
value of d(^/b)/dt is determined, the kinetic energy of the plate will decrease as ?(t)/b 
decreases and the time rates of curvature k r and kg will be nonnegative also. Hence, 
solving equation (C4) for the largest value of ?(0)/b for which d(£/b)/dt is negative 
yields 



For values of ^ > 1, the hinge circle remains at the outer edge of the sheet since the 
simple support boundary condition M r (b,t) = 0 (eq. (C5)) prohibits regime A from 
encompassing the edge of the plate and passing beyond. Hence for ^ g 1, 


m 

b 


= 1 


(C8) 


The value of t = tj, at the time the hinge circle vanishes at the center of the sheet, 
can be determined from equation (C6) by first evaluating Cj ^by setting t = 0 and 
k = ( e( l s ' (C7) or (C8))) and then setting £ = 0 to yield 


1 

! 

I 
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24M ° rtl = -2 + afef + ,1 - (A 


MV 0 b 


2 l l 


2 - 6 (t ) 2 + (t ) 4 


0 g s < l] 


and 


24M r 


MV 0 b 


2 t 1 = -2 + 8(|) +,/ 1 -[£ 


12 - 5 I 


3f£T sin -1 | 

a 
b 


3(£ ) sin"l — 
\b / a 


(C9) 


#>1 


(CIO) 


In general, for ^ = 4, equation (C6) becomes 


24M, 


/iV 0 b2 


%t(4) = 


1 “ (iT7b) 


ia\ 2 9 /a\ 2 

-2 + 8fgJ - 2| + 2| 2 _ 3y±j | 

/ a \ 3 i / A \ /a\ 2 24M 0 

+ 3(£ sin _1 (-i-] + 2 - 8(f + 

v/ \a/b/ v/ uV n b 2 


(Cll) 


The deflection of the plate with a hinge circle of finite radius can now be deter- 
mined since the velocities of deformation have been established. The deflection of the 

NO) 

surrounding annulus of the plastically deforming sheet, from — — — S 4 = 1, which has 
always remained in plastic regime AB (eq. (C2)) is 


w 



Vr 


!-iy 2 r 


1/2 




>-§ 


dt 



(Cl 2) 


The central portion of the plastic sheet contains sheet elements that are first located in 
the plastic regime A and then, as the hinge -circle passes over them, in the plastic 
regime AB. The deflection is, therefore, dependent upon the time integration of two 
velocities: one velocity when the point in question is in regime A (eq. (Cl)), and the 
second velocity (after time t(4) when the hinge circle passes over that point (£ = 
when the point is in regime AB (eq. (C2)). Hence, ' 


w(!,ti) = vJ 1 -( 1 ) £ 2 t(*) 




(Cl 3) 
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Expressing dt in terms of d(£/b) 
yields 
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(eq. (C4)) and integrating equations (C12) and (Cl 3) 


w 


M = 


M V 0 2 b 2 

24M n 


(1 - €) 


f ) 2 + (h 


(C14) 


w 


(«.‘l) - v o/l -( + ^£-<1 - «) ai + 3{2 - (§ J { 2 - 5 3 


where t(|) is given by equation (Cll). 




?( 0 )' 


(Cl 5) 


Solution With Hinge Circle of Zero Radius 

When the hinge circle radius shrinks to zero at t = tj and the entire plate is in the 
plastic regime AB, the previous analysis must be terminated because the acceleration of 
the plate vanishes and the equation of motion (see eq. (C3)) reduces to 

^ (« M r) = Mo 
or 

M r = M 0 + ^ (Cl 6) 

where C3 is the constant of integration. Since M r = M 0 at ij = 0, C3 = 0 and M r 
remains at the constant yield-moment value of M 0 in contradiction to the prescribed 
plasticity condition in regime AB (ref. 4) and the simply supported edge condition at | = 1 
(eq. (C5)). It can be shown that the following expression for sheet deflection satisfies the 
plasticity and boundary conditions of the sheet, once the hinge-circle radius has shrunk to 
zero: 


w(£,t) = (p (t) (1 - 0 + w(|,ti) (t ^ tj) (Cl 7) 

where w(|,ti) is the deflection at t = ti, given by equations (C14) and (Cl 5), and <p( t) 
is an arbitrary function of time. 

Again, the acceleration is obtained by differentiation of equation (Cl 7) and is substi- 
tuted into equation (C3) with Mg = M 0 . The subsequent integration of equation (C3) and 
application of the boundary conditions: 
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yield 
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M r = M 0 
M r = 0 


({ = 0)1 

> (Cl 8) 

({ = D 

-s 



(Cl 9) 


which, upon integration by utilizing the conditions evaluated at t = tj, 


becomes 



(C20) 


■»<*> = - 5(‘" tl ) 2 + v °( t "' 1 ) 


(C21) 


Hence, the deflection for t = tj is 


w(Z,t) = + (1 - I) - t : ) 2 + V 0 (t - tjV (t 2 tj) (C22) 

- _| 


The time t = t* at which the entire plate comes to rest can be found by differentiating 
equation (C22) to obtain the velocity, by equating the resulting expression for velocity to 
zero, and then solving for t*. This procedure yields 


12M 0 t% 

MV 0 b2 


l2M n 



MV 0 b 2 


(C23) 


and the permanent deformation shape, denoted as w*(£,t*) for t = t*, can then be 
written as 


w*({,t*) = w^tj) + (1 


- 0 


MV 0 2 b 2 
24M 0 


(C24) 
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The permanent center deformation becomes (eqs. (C15) and (C24) with | = 0) 


w*(0,t*) 


6 = 


jHV 0 2 b 2 

24M 0 




(C25) 
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TABLE I.- TARGET AND PROJECTILE MATERIAL PROPERTIES 
AND VALUES OF THE MOMENTUM CONSTANT \p 


(a) Target and projectile material properties 


Property 

U.S. unit 

SI unit 

Target material (2024-T4 aluminum) 


Yield stress, a Q 

68 000 psi 

0.47 GN/m2 

Modulus of elasticity, E . . . 

10.6 x 10® psi 

73.1 GN/m2 

Density, pj. 

5.37 slugs /ft ^ 

2767 kg/m3 

O 

II 

r+ 

16 900 ft/sec 

5.15 km/s 

Projectile 


Density, p p 

1.78 slugs /ft 3 

917 kg/m 3 

Specific gravity 

0.916 

0.916 


Target 

thickness, 

h 


in. 

0.063 

.04 


cm 

0.160 

.10 


(b) Values of the momentum constant i// 
computed from equation (7) 

Values of \p at various projectile radii 


a = 0.50 in. (1.27 cm) 
0.70 
.69 


a = 0.25 in. (0.64 cm) 
0.71 
.70 


a = 0.125 in. (0.32 cm) 


0.72 

.71 
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Nondimensional 

moment 



Figure 3.- Typical moment distributions in central portion of infinite 



plate, t = 1.5. 


Dimensionless radius 
plastic region , 7 

d 



Figure 4.- Radius of plastic region as a function of projectile velocity. 




100 1000 10000 


Projectile velocity, g , ft/s 

Center deflection 6 as a function of projectile velocity g 0 . 0.063-inch-thick (0.160-cm) 2024-T4 aluminum sheet; 6 





(a) Impact velocity, 330 ft/sec (0.10 km/s). 



(b) Impact velocity, 500 ft/sec (0.15 km/s). 



(c) Impact velocity, 1100 ft/sec (0.34 km/s). 

Figure 6.- Variation in permanent deflection of 2024-T4 aluminum sheet impacted with 1-inch-diameter projectile 

at various velocities, a = 0.50 in. (1.27 cm). 


Center def lectior 
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Figure 7.- Center deflection 6 as a function of projectile velocity g 0 . 0.04-inch-thick (0.10-cm) 2024-T4 aluminum sheet; 6 ■= g 0 
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